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Abstract 

In this article, we calculate the contributions of the vacuum condensates up to dimension-6 
including the (D{a a ) corrections to the quark condensates in the operator product expansion, 
then study the masses and decay constants of the pseudoscalar, scalar, vector and axial- 
vector heavy-light mesons with the QCD sum rules in a systematic way. The masses of the 
observed mesons (£>,£>*), (D S ,D* S ), (I>3(2400), £>i(2430)), (D* 0 (2317), £> s i(2460)), 

( B S ,B *) can be well reproduced, while the predictions for the masses of the (Bg,Bi) and 
(B* 0 ,B s i) can be confronted with the experimental data in the future. We obtain the decay 
constants of the pseudoscalar, scalar, vector and axial-vector heavy-light mesons, which have 
many phenomenological applications in studying the semi-leptonic and leptonic decays of the 
heavy-light mesons. 

PACS number: 13.20.Fc, 13.20.He 
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1 Introduction 


The charged heavy-light mesons can decay to a charged lepton pair £ + ve through a virtual W + 
boson. Those leptonic decays are excellent subjects in studying the CKM matrix elements and 
serve as a powerful probe of new physics beyond the standard model in a complementary way to the 
direct searches. For example, the decay widths of the pseudoscalar (P) and vector (V) heavy-light 
mesons can be written as 


r(P -A tv) = |/Xm P (l-|) \V qiq2 \ 2 , 

T{V^lv) = !-^r) ( 1 + ^)l^! 2 


(1) 


in the lowest order approximation, where the mp/y and fp/v are the masses and decay constants, 
respectively, the me is the £ mass, the V qiq2 is the CKM matrix element between the constituent 
quarks q±q 2 , and the Gf is the Fermi coupling constant. If we take the CKM matrix element V qiq2 
and the branching fractions of the leptonic decays from the CLEO, BaBar, Belle collaborations as 
input parameters, then the average values fp = (204.6 ± 5.0) MeV, fp a = (257.5 ± 4.6) MeV and 
fn s /fn = 1.258 ± 0.038 are obtained pQ. It is difficult to reproduce the three values consistently 
in theoretical calculations, such as the QCD sum rules mmmm and lattice QCD [6j[7l|8]. The 
discrepancies between the theoretical values and experimental data maybe signal some new physics 
beyond the standard model [5j. In Ref. M, we observe that if we take into account the O(a^) 
corrections to the perturbative terms and the 0(a s ) corrections to the quark condensate terms 
and choose the pole masses, the predictions fp = (211 ± 14) MeV, fp s = (258 ± 13) MeV and 
Id s /Jd = 1-22 ± 0.08 are in excellent agreement with the experimental data [T. 

In the QCD sum rules for the heavy-light mesons, the Wilson coefficients of the vacuum con¬ 
densates at the operator product expansion side from different references differ from each other in 
one way or the other according to the different approximations [2 0 [nil HUE]. In this article, we 
recalculate the contributions of the vacuum condensates up to dimension-6, including the one-loop 
corrections to the quark condensates, and take into account the terms neglected in previous works, 

1 E-mail,zgwang@aliyun.com. 


1 




then study the masses and decay constants of the pseudoscalar, scalar, vector and axial-vector 
heavy-light mesons in a systematic way. 

There have been many theoretical works on the decay constants of the heavy-light mesons, 
such as the QCD sum rules 0 02 m 0 0 US HU H2 Q21 Ull [HI Q21H3 HH1HH1 EHE], the lattice 
QCD [6] [7] [8] E2 EH EH EH EH EZ] ; the Bethe-Salpeter equation [28l [29] . the relativistic potential 
model mm, the field-correlator method [32] , the light-front quark model [3H |31] , the chiral 
extrapolation [35], the extended chiral-quark model [36], the constituent quark model m , etc. 

The article is arranged as follows: we derive the QCD sum rules for the masses and decay 
constants of the heavy-light mesons in Sect.2; in Sect.3, we present the numerical results and 
discussions; and Sect.4 is reserved for our conclusions. 


2 QCD sum rules for the heavy-light mesons 

In the following, we write down the two-point correlation functions II0/5 (p) and H^j A {p) in the 
QCD sum rules, 


n0/5 (p) 

= i J d 4 xe ipx {0\T { J 0/5 (a:)4 /5 (0)} |0), 

(2) 

n V}a(p) 

= i j d 4 xe ipx (0\T{j» /A (x)J"} A ( 0)} 0), 

(3) 

Jo( x) 

= Q(x)q(x), 


J 5 (x) 

= Q(x)i'j5q(x) , 


Jv(x) 

= Q(.x)Yq{x) , 


Ja(x) 

= Qixh^^qix), 

(4) 


where the currents J${x), Jq(x), Jy(x) and Jfix) interpolate the pseudoscalar, scalar, vector and 
axial-vector heavy-light mesons, respectively, Q = c,b and q = u,d, s. We can insert a complete set 
of intermediate hadronic states with the same quantum numbers as the current operators ^(x), 
Jo{x), Jy{x) and J A (x) into the correlation functions n 0 / 5 (p) and II ^f/ A (p) to obtain the hadronic 
representation [38][39]. After isolating the ground state contributions from the pseudoscalar, scalar, 
vector and axial-vector heavy-light mesons, we get the following results, 


n 0 (p) 

n 5 (p) 

K/a(p) 


fs m s 

9 9 

771 g ~ P 


fp m 


('rriQ + m q ) 2 (rn 2 p - p 2 ) 

f2 

J i 


IV/A m V/A 


m v/ A - P 


+ 


p^p v 


pPp" 


n v/A(p)[-gpu+ p2 


where the decay constants fs/p/v/A are defined by 


<0|J 0 (0)|S(p)) 

<0|J 5 (0)|P(p)> 

(0\J» /A mV/A(p)) 


= fsms, 
fpm 2 P 
mQ + m q ’ 

= fv/Amy/A^ , 


(5) 

( 6 ) 

(7) 


( 8 ) 


the e p are the polarization vectors of the vector and axial-vector mesons. 
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Figure 1: The diagram contributes to the quark condensate (qq). 



Figure 2: The perturbative 0(a s ) corrections to the quark condensate (qq). 


Now we carry out the operator product expansion at large space-like region P 2 = —p 2 . The 
analytical expressions of the perturbative 0(a s ) corrections to the perturbative terms for all the 
correlation functions mm and the semi-analytical expressions of the perturbative 0(a 2 ) correc¬ 
tions to the perturbative terms for the pseudoscalar current’s correlation functions (19] are available 
now. We take into account those analytical and semi-analytical expressions directly nzmamsi; 
and recalculate the contributions of the vacuum condensates, i.e. we calculate the Feynman dia¬ 
grams shown in Figs. 1-5, where the solid and dashed lines denote the light and heavy quark lines, 
respectively, the wave line denotes the gluon line. In calculating the diagrams in Fig.2, we correct 
the minor errors in Ref. uni, where the quark condensate -^jrr in the full light-quark propagators 
is replaced with , the D is the dimension of the space-time. A minor error occurs when there 

exist divergences, such a step should be deleted, i.e. the quark condensate survives in the 
.D-dimension. In Ref. j40j. we correct the minor errors and improve the calculations, and obtain 
the correct expressions. Furthermore, we obtain the perturbative 0(a s ) corrections to the quark 
condensate terms for the vector and axial-vector currents. 



Figure 3: The diagrams contribute to the mixed condensate ( qg s aGq). 
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Figure 4: The diagrams contribute to the gluon condensate ( as ^ G ) and three-gluon condensate 
(. 9sGGG ). 
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Once analytical expressions of the QCD spectral densities are obtained, then we can take 
the quark-hadron duality below the continuum thresholds and perform the Borel transforms with 
respect to the variable P 2 = — p 2 to obtain the QCD sum rules, 


where 


fp m p c::p 

(m Q + m. q y \ T 2 ) 
/im|exp(-^|) 

fv m v exp 

f\ m \ exp 


Bt II5 , 


BtBq , 


Bt By , 


BtBa , 


B t n 5 = B t U ° 5 + B t U 3 5 + Ij T 115 + /irIJ 5 + B T nf , 

BtBq = BtB^ rriQ^—mQ •> 

BtBv — BtB^ + Bt^y + Bt^y -)- + H^IIy , 

= HtHv ImQ—t—mQ i 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


(13) 

(14) 


B t n° 
B t U 3 5 


B t Ilj 
Bt n« 

s T nj? 
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Bt K = 


Bt 111 1 - = 


BtK 


B t H 5 v = 


B t n 6 


v 



m %{ ( l9svGq) m q (qg s aGq) 


207T a s (qq) 2 
81 T 2 


m Q 

' 768T 4 


1 - 


1 + 


'"Q 

rp2 


Q 


Q 


T 2 5 T 4 
55 mg 


exp - 


m Q\ (glGGG) 
T 2 n 2 


1 


1152T 2 


1152mQ 


+ 


1 


m r 


4608T 4 192T 2 


, , l Q m Q \, rn 2 Q g 2 

1+^2 27 1 4 j lo § y 4 


exp — 


^2 


(24) 


Rv{x) 


Q 3 1 _ X 1 _ X X 

- + 2Li 2 (x) + logo; log(l — x) — - log-log(l - x) + x log---logx, 

4 lx x 1 — x 


+ 2Li 2 (a;) + logxlog(l - x) - ^ log-—- - log(l - x) + Z log-—- - —— 
4 lx x 1 — x 


(3 + a;)(l — x) . 1 — x 

■log 


2x 


2 T x 


x (2 + x){\ — x) 2 


log a; — 


2a; 


5 + 2a; 

2 + x (2 + a;)(l — x) 


logs 

(25) 


r°° i 

T(0,x) = e x dt -e 1 , 

Jo t + x 

Li 2 (a;) = - [ dG- log(l - t) , (26) 

Jo t 

and the so are the continuum threshold parameters. The perturbative 0(a s ) corrections i?s(a;) and 
Rv[x) are taken from Refs.T7 s 18|. We can also take into account the semi-analytical perturbative 
0(a 2 ) corrections to the perturbative terms for the BtB° , 

i (' v) 2 / l ds {T R2sFFM + 4R2sFA M + 2 ~y R2sFL M +\ R2sFH M } ex P (-^2 ) > (27) 

where the R2sFF[u], R2sFA[u], R2sFL[u] and R2sFH[r>] with the variable v = ^1 — / ^1 + 

are mathematical functions defined at the energy-scale of the pole mass g = m c , here the ni counts 
the number of massless quarks [19] . 

We can derive Eqs.(9-12) with respect to 1/T 2 , then eliminate the decay constants fs/p/v/A 
to obtain the QCD sum rules for the masses. 


_ 2 _ d(l/T' z ) BtIIo/ 5 /V/A 

B T U 0/5/v/A 


(28) 
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Once the masses ms/p/v/A are obtained, we can take them as input parameters and obtain the 
decay constants from the QCD sum rules in Eqs.(9-12). 

In the case of the light-quark currents, the perturbative 0(a s ) corrections to the perturbative 
terms amount to multiplying the factors 1 ± qp' ~ 1 + 3.67-^r and 1 + ^ to the perturbative 
terms in the correlation functions for the pseudoscalar (scalar) and vector (axial-vector) currents, 
respectively [35]. In the present case, if we take the approximation /i 2 = to 2 = T 2 , the perturbative 
0{a s ) corrections to the quark condensate terms amount to multiplying the factors 1 ±3.47^ and 
1 ± 0.94-^f- to the quark condensate terms in the correlation functions for the pseudoscalar (scalar) 
and vector (axial-vector) currents, respectively. The analogous 0(a s ) corrections indicate that the 
present calculations are reliable. 


3 Numerical results and discussions 


In the heavy quark limit, the heavy-light mesons Qq can be classified in doublets according to the to¬ 
tal angular momentum of the light antiquark Sf . S£ = Sq+L , where the Sq and L are the spin and or¬ 
bital angular momentum of the light antiquark, respectively. The spin doublets ( D , D*), ( D s , D*), 
(D* (2400),Di(2430)), (D* 0 (2317),D s i(2460)), (B S ,B* S ) have been observed, the masses 

are m D ± = (1869.5 ± 0.4) MeV, m D o = (1864.84 ±0.07) MeV, m D . (2010) ± = (2010.26 ±0.07) MeV, 
ttid* ( 2007 )° ~ (2006.96 ± 0.10) MeV, r77.£)*f24oo± = (2318 ± 29) MeV, tyi b *(24004 = (2403 ± 14 ± 
35) MeV") m Dl ( 2430) o = (2427± 26 ± 25) MeV= (1969.0± 1.4) kv) m"= (2112.1 ± 
0.4)MeV, m D . o(23 i 7 )± = (2318.0±1.0)MeV, m Da ’ {2i60) ± = (2459.6±0.9)MeV, m B ± = (5279.25± 
0.26) MeV, m B o = (5279.55 ± 0.26) MeV, rn B - = (5325.2 ± 0.4) MeV, m Be = (5366.7 ± 0.4) MeV, 
m B » = (5415.8 ± 1.5) MeV from the Particle Data Group pQ. The spin doublets (Z?g,i?i) and 
(B* 0 , B s i) have not been observed yet. The doublet (D(2550), D(2600)) or (Dj(2580), D}(2650)) is 
tentatively identified as the first radial excited state of the doublet ( D , D*), the doublet (?, D* sl (2700)) 
is tentatively identified as the first radial excited state of the doublet (D S ,D*( 2112)) [44] . 

We take the values ^/sq = m gr ± (0.4 — 0.8) GeV as guides, here the gr denotes the ground 
states, and search for the optimal threshold parameters so and Borel parameters T 2 to satisfy the 
following criteria: 

• Pole dominance at the phenomenological side; 

• Convergence of the operator product expansion; 

• Appearance of the Borel platforms; 

• Reappearance of experimental values of the ground state heavy meson masses. 

The contributions of the ground states can be fully taken into account by choosing the threshold 
parameters yGo = m gr ± (0.4 — 0.8) GeV. The contaminations of the excited states are very small 
if there are some contaminations, we expect that the couplings of the currents to the excited states 
are more weak than that to the ground states. For example, the decay constants of the pseudoscalar 
mesons 7r(140) and 7r(1800) have the hierarchy /^moo) <C /„-(i 4 o) from the lattice QCD [45], the 
QCD sum rules [43], or from the experimental data [44] . 

The vacuum condensates are taken to be the standard values (uu) = (dd) = — (0.24±0.01 GeV) 3 , 
(ss) = (0.8 ± 0.1)(uu), ( qg s aGq) = m^{qq}, mg = (0.8 ± 0.1) GeV 2 , { as ^ G ) = (0.33 GeV) 4 , 
{g^GGG) = 0.045 GeV 6 at the energy scale jj. = 1 GeV J38] iSSV The quark condensates and mixed 


quark condensates evolve with the renormalization group equation, (qq)(n) 


(qq)(Q) 


MQ) 1 9 
«s(/d 


and ( qg s aGq)(ii ) = ( qg s aGq)(Q ) 

In the article, we take the MS masses mb{rrib ) = (4.18 ± 0.03) GeV, m c {m c ) = (1.275 ± 
0.025) GeV and m s (g = 2 GeV) = (0.095 ± 0.005) GeV from the Particle Data Group [I], and 
take into account the energy-scale dependence of the MS masses from the renormalization group 
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equation, 


rrib(^) 

m c (n) 

m s {n) 

m u / d (n) 

a s (/x) 


rrib{mb) 

m c (m c ) 


a 


(m) 


a s (m b ) 

a s (/u) 


m s (2GeV) 
«V d (lGeV) 
1 - 


_a s (m c )_ 

u s {n) 5 


_a s (2GeV)J ’ 

Mm) 15 


1 


a s (lGeV)J ’ 

&ilogt & 2 (log 2 f - logt - 1) + b Q b 2 


b 2 0 t 


bfc 2 


(29) 


where t = log ^ 
and 339 MeV for the flavors 


7 33—2rif 7 153—19 nt i 

b 0 = 127T ' b 1 = 24^ V & 2 = 

n/ = 5, 4 and 3, respectively [F. 


2S57— 5033 r? 325 „2 

—— 1 28 X 27 7 , a = 213 MeV, 296 MeV 
Furthermore, we obtain the values 


m u = rrid = 6 MeV from the Gell-Mann-Oakes-Renner relation at the energy scale /j, = 1 GeV. 

In this article, we choose the MS masses by setting m = and take the perturbative 

0(a s ) corrections to the perturbative terms. In other words, we take the R$/v () only. In 


calculations, we take nj = 3 and /j, d/d * = s/ m|, — m 2 ~ 1 GeV for the S-wave mesons D and D*; 
rif = 4 and ^b/b* = \Jm 2 B ~ m b ~ 2.5 GeV for the S-wave mesons B and B*. If we count the 
contribution of the additional P-wave as 0.5 GeV, then £ 1 d*/d 1 = 1-5 GeV and ^b^/Bx = 3.0 GeV. 
On the other hand, we take into account the SU{3) breaking effect, which is supposed to be 
100 MeV for the light quarks, then n d 3 /d* =1-1 GeV, Hb„/B‘ =2.6 GeV, Hd* 0 /d b1 =1-6 GeV and 
Hb* 0 /b s1 = 3.1 GeV. Those energy scales work well. 

The continuum threshold parameters, Borel parameters, pole contributions are shown explicitly 
in Table 1. From Table 1, we can see that the pole dominance can be satisfied. On the other hand, 
the dominant contributions come from the perturbative terms and the quark condensate terms, so 
we expect to obtain reliable predictions. 

After taking into account the uncertainties of the input parameters, we obtain the values of the 
masses and decay constants of the heavy-light mesons, which are shown in Figs.6-9 and Table 1. 
From the figures, we can see that the masses and decays constants are rather stable with variations 
of the Boral parameters T 2 , the predictions are reasonable. 

From Table 1, we can see that the experimental values of the masses of the observed mesons 


(D,D*), (D S ,D*), (I>5(2400), £>4(2430)), (H s * 0 (2317), £> sl (2460)), (B S ,B*) can be well 

reproduced. The masses of the (Bq,B i) and (B* 0 ,B S i) vary in rather large ranges from different 
theoretical approaches, ms* = (5.53 — 5.76) GeV, ms ± = (5.58 — 5.78) GeV, 77 is* o = (5.63 — 
5.83) GeV, niB Bl = (5.67 — 5.86) GeV, for a comprehensive review, one can consult Ref. [451. The 
present predictions to_bj = (5.72±0.05) GeV, mB 1 = (5.74±0.05) GeV, t?xb* o = (5.70±0.06) GeV, 
w_b s1 = (5.76 ± 0.06) GeV are compatible with those values. 

The thresholds of the DK, D*K, BK and B*K states are vidk = 2.36 GeV, mn^K = 
2.50 GeV, nr bk = 5.78 GeV and tub^k = 5.82 GeV, respectively. The £>* 0 (2317) and £> s i(2460) he 
below the thresholds vtidk and respectively, the Okubo-Zweig-Iizuka allowed strong decays 

D* 0 (2317) —> DK and £> s i(2460) —> D*K are kinematically forbidden, the widths of the £>* 0 (2317) 
and £> s i(2460) are very narrow. According to the present predictions mB* 0 = (5.70 ± 0.06) GeV 
and niB Bl = (5.76 ± 0.06) GeV, the B* 0 and B si also lie below the corresponding BI\ and 
B*I\ thresholds, respectively. The strong decays B* 0 —> BK and B s i —»• B*K are kinemati¬ 
cally forbidden, the P-wave heavy mesons B* 0 and B s \ can decay through the isospin violation 
precesses B* 0 —► B s r] —> B s n° and B s i —» B*rj B*tt° respectively or through the radiative 
decays [46]. The 77 and 7 r° transition matrix is very small according to Dashen’s theorem [47], 
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T 2 , the A, B, C , D , E, F, G and H denote the mesons D , D s , D*, D*, Dq, D* 0l D\ and D s i, 
respectively. 
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Figure 9: The decay constants of the bottom mesons with variations of the Borel parameters 
T 2 , the A, B , C, D , E , F, G and H denote the mesons B , B s , B*, B*, Bq, B * 0 , B\ and B s \, 
respectively. 
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T 2 ( GeV 2 ) 

s 0 (GeV 2 ) 

pole 

ms/p/v/A{GeV) 

fs/p/v/A{ MeV) 

D 

1.2- 1.8 

6.2 ±0.5 

(67 - 93)% 

1.87 ±0.10 

208 ± 10 

D s 

1.3- 1.9 

7.3 ±0.5 

(77 - 96)% 

1.97 ± 0.10 

240 ± 10 

D* 

1.9-2.5 

6.4 ±0.5 

(51 - 76)% 

2.01 ±0.08 

263 ±21 

D* 

1.8- 2.4 

7.5 ±0.5 

(66 - 87)% 

2.11 ±0.07 

308 ±21 

D* 

2.4- 3.0 

8.3 ±0.5 

(59 - 78)% 

2.40 ±0.05 

373 ± 19 

D* s o 

2.1 - 2.7 

7.4 ±0.5 

(55 - 77)% 

2.32 ±0.05 

333 ± 20 

D\ 

2.9- 3.5 

8.6 ±0.5 

(52 - 70)% 

2.42 ±0.05 

332 ± 18 

D s i 

2.7- 3.3 

9.3 ±0.5 

(61 - 78)% 

2.46 ±0.06 

245 ± 17 

B 

5.5-6.5 

34.0 ± 1.0 

(44 - 63)% 

5.28 ±0.07 

194 ±15 

B s 

5.3-6.3 

36.0 ± 1.0 

(56 - 74)% 

5.37 ±0.07 

231 ±16 

B* 

6.0- 7.0 

34.5 ± 1.0 

(45 - 62)% 

5.32 ±0.06 

213 ±18 

B* s 

6.1 - 7.1 

36.5 ± 1.0 

(53 - 69)% 

5.42 ±0.06 

255 ± 19 

B*o 

6.3- 7.3 

40.0 ± 1.0 

(60 - 75)% 

5.72 ±0.05 

281 ±14 

R* 

^sO 

6.5 - 7.5 

40.0 ± 1.0 

(60 - 74)% 

5.70 ±0.06 

274 ± 13 

Hi 

6.7- 7.7 

41.0 ± 1.0 

(62 - 75)% 

5.74 ±0.05 

335 ± 18 

B s i 

7.1 - 8.1 

42.0 ± 1.0 

(63 - 76)% 

5.76 ±0.06 

348 ± 18 


Table 1: The Borel parameters, continuum threshold parameters, pole contributions, masses and 
decay constants of the heavy-light mesons. 

t V ir = (7r°|H|r?) = —0.003 GeV 2 , the P-wave bottomed mesons B* 0 and B s i, just like their charmed 
cousins D* 0 (2317) and D s i(2460), maybe very narrow [IS]. The present predictions are consistent 
with our previous work [13], but the analysis is refined by including more terms in the operator 
product expansion. 

The values of the decay constants of the pseudoscalar mesons are slightly different from the ones 
in our previous work [101 . In Table 2, we compare the present predictions to the experimental data 
and other theoretical calculations, such the QCD sum rules (QCDSR) [2:, 3] [I, ;5J jTIj and lattice 
QCD (LQCD) [U El [8]. The present predictions fn = (208 ± 10) MeV and fs = (194 ± 15) MeV 
are consistent with the experimental data within uncertainties, while the prediction fjj s = (240 ± 
10) MeV is lies below the lower bound of the experimental value fo e = (257.5 ± 4.6) MeV pQ. 
We take the MS mass m c (n) and truncate the perturbative corrections to the order 0(a s ), the 
experimental values of the fn, /d„ and fnj Id cannot be reproduced consistently by the QCD sum 
rules. The existence of a charged Higgs boson or any other charged object beyond the standard 
model would modify the decay rates, see Eq.(l), therefore modify the values of the decay constants, 
for example, the leptonic decay widths are modified in two-Higgs-doublet models @5]. If the 
predictions of the f D , f Ds and /d„//d based on the QCD sum rules are close to the true values, 
new physics beyond the standard model are favored so as to smear the discrepancies between the 
theoretical calculations and experimental data. 

The analytical expression of the perturbative 0{a s ) corrections R 5 is well known [T7j, 

while the semi-analytical perturbative 0(a 2 ) corrections are presented as mathematical functions 

R2sFF[z>], R2sFA[w], R2sFL[ri] and R2sFH[ri] with the variable v = ^1 — / (l + at the 

energy-scale of the heavy quark pole mass /r = toq ll9j . The analytical expressions of the terms 

2 2 

which contain logarithms such as log -±j-, log cannot be recovered, it is unreasonable to take 
other energy scale besides toq. Now we choose the pole masses ?tiq and take into account the 
semi-analytical perturbative 0(a 2 ) corrections by setting rif = 4 and fi = m c for the D ( D s ) 
meson and rif = 5 and fx = mb for the B ( B s ) meson. 

The on-shell quark propagators have no infrared divergences in perturbation theory, which 
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provides a perturbative definition of the quark masses. The full quark propagators have no poles 
because the quarks are confined, so the pole masses cannot be defined outside of perturbation 
theory. Furthermore, the pole masses cannot be used to arbitrarily high accuracy because non- 
perturbative infrared effects in QCD. We choose the pole masses just because the semi-analytical 
perturbative O(a^) corrections are calculated by taking the pole mass rriQ and setting the energy 
scale to be /i = toq IS]. The contributions of the u, d masses are tiny and can be neglected 
safely. In calculations, we set the pole masses m u = rrid = 0, m s = 150 MeV, and observe that the 
masses of the heavy pseudoscalar mesons increase monotonously with increase of the pole masses, 
the values of the pole masses m c = 1.44 GeV and mb = 4.67 GeV can lead to satisfactory values 
by choosing reasonable Borel parameters and threshold parameters. Those pole masses are differ¬ 
ent from the MS masses, for example, m c (/z = 1 GeV) = 1.39 GeV, mb(n = 1 GeV) = 6.07 GeV, 
m c (n = 2 GeV) = 1.13 GeV, mt(/i = 2 GeV) = 4.87 GeV from Eq.(29). The pole masses are energy 
scale independent, therefore the energy scale dependence of the QCD spectral densities originate 
only from the vacuum condensates. 

The Borel parameters, continuum threshold parameters, pole contributions, and the resulting 
masses and decay constants of the heavy pseudoscalar mesons are shown in Table 3, the values are 
slightly different from the ones in our previous work m • From Table 1 and Table 3, we can see 
that the present predictions fp> = (210± 11) MeV, fp s = (259± 10) MeV and fg = (192± 13) MeV 
are in excellent agreement with the experimental data within uncertainties 0 • The ratio fnjfn = 
1.23 ± 0.07 is also in excellent agreement with the experimental data Jd s /Id = 1.258 ± 0.038 if], 
which indicates that the perturbative O(otg) corrections should be taken into account. However, 
the pole masses toq and energy scales n = mQ have be chosen, as the semi-analytical expressions 
are obtained at such conditions. In this case, new physics beyond the standard model are not 
favored, as the agreements between the experimental data and present theoretical calculations are 
already excellent. 

In Table 4, we compare the present predictions for the decay constants of the heavy vector 
mesons to other theoretical calculations, such as the QCD sum rules [3 EE GU, lattice QCD 
[23 [231 M HE the relativistic potential model (RPM) [30], the field-correlator method (FCM) 
[32] . and the light-front quark model [34]. From the table, we can see that the predictions differ 
from each other in one way or the other. In Table 5, we compare the present predictions for the 
decay constants of the heavy scalar mesons to the ones from the QCD sum rules m and lattice 
QCD [25] , From the table, we can see that the predictions are consistent with the ones from lattice 
calculations but differ greatly from the ones from the QCD sum rules. 

If we turn off the perturbative 0(a s ) corrections to the quark condensates and choose the same 
parameters, such as the MS masses, Borel parameters and continuum threshold parameters, etc, 
the masses and decay constants undergo reduction or increment in a definite way according to the 
spin and parity, see Table 6. From the table, we can see that the mass-shifts of the D-mesons 
with J p = 0^ are larger than 40 MeV, while the shifts of the masses and decay constants of all 
the H-mesons are small and can be neglected. We can re-choose the Borel windows to warrant the 
mass-shifts Sms/p/v/A = 0, and account for the net effects by the shifts of the decay constants 
fifs/p/v/Ai which are shown the bracket in Table 6. From the table, we can see that the largest 
shift of the decay constant Sfp> = —11 MeV, which exceeds the total uncertainty of the decay 
constant 5fp = ±10 MeV (see Table 1), the shifts of the decay constants of the D-mesons with 
J p = O^, 1 — are larger than 5 MeV, while for other mesons, the shifts of the decay constants 
\5f\ < 4 MeV. All in all, we should take into account the perturbative 0(a s ) corrections to the 
quark condensates in a comprehensive study. 


4 Conclusion 

In this article, we calculate the contributions of the vacuum condensates up to dinrension-6, in 
including the 0(a s ) corrections to the quark condensates, in the operator product expansion. Then 
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}d{ MeV) 

Id s (MeV) 

M MeV) 

fsA MeV) 

/dJ Id 

fBs/ 7b 

Expt pQ 

204.6 ± 5.0 

257.5 ±4.6 

190.6 ±4.7 


1.258 ±0.038 


QCDSR (2j 

177 ±21 

205 ± 22 

178 ±14 

200 ± 14 

1.16 ±0.16 

1.12 ± 0.11 

QCDSR 0 

204 ±6 

246 ±6 

207 ±8 

234 ±5 

1.21 ±0.04 

1.14 ±0.03 

QCDSR g] 

206.2 ± 7.3 

245.3 ± 15.7 

193.4 ± 12.3 

232.5 ± 18.6 

1.193 ±0.025 

1.203 ±0.020 

QCDSR 0 

20 F^ 

9 QQ + 13 

^00 23 

207 +iY 

zu '-09 

2421% 

1 1 Q+U-04 

l-lo o.o5 

1 1 7+U.U3 

1 -0.04 

QCDSR [g 



186 ± 14 

222 ± 12 


1.19 ±0.09 

LQCD |5j 

197 ±9 

244 ±8 



1.24 ±0.03 


LQCD [7] 

213 ±4 

248.0 ±2.5 

191 ± 9 

228 ± 10 

1.164 ±0.018 

1.188 ±0.018 

LQCD Ej 

218.9 ± 11.3 

260.1 ± 10.8 

196.9 ±8.9 

242.0 ±9.5 

1.188 ±0.025 

1.229 ±0.026 

This work 

208 ± 10 

240 ± 10 

194 ± 15 

231 ±16 

1.15 ±0.06 

1.19 ±0.10 

This work* 

210 ± 11 

259 ± 10 

192 ± 13 

230 ± 13 

1.23 ±0.07 

1.20 ±0.09 


Table 2: The decay constants of the heavy pseudoscalar mesons from the experimental data, the 
QCD sum rules and lattice QCD, the superscript star * denotes that the pole masses are chosen 
and perturbative 0 ( 0 ^) corrections are taken into account. 



T 2 ( GeVQ 

so(GeV") 

pole 

mp( GeV) 

fp( MeV) 

D 

1.4- 2.0 

5.5 ±0.5 

(55 - 85)% 

1.87 ±0.06 

210 ± 11 

D s 

1.0- 1.6 

7.4 ±0.5 

(86 - 98)% 

1.97 ±0.07 

259 ± 10 

B 

4.1-4.9 

33.0 ± 1.0 

(55 - 75)% 

5.28 ±0.04 

192 ± 13 

B s 

4.4- 5.2 

35.0 ± 1.0 

(61 - 79)% 

5.37 ±0.04 

230 ± 13 


Table 3: The Borel parameters, continuum threshold parameters, pole contributions, masses and 
decay constants of the heavy pseudoscalar mesons when the perturbative O(cts) corrections are 
taken into account. 



f D -{ MeV) 

/dj (MeV) 

/s*(MeV) 

/ B: (MeV) 

QCDSR [5] 

242^12 

293+]^ 

210™ 

25l^ 

qcdsr pn 

252.2 ±22.3 ±4 

305.5 ±26.8 ±5 

181.8 ±13.1 ±4 

213.6 ± 18.2 ±6 

QCDSR [21] 

250 ± 11 

270 ± 19 

209 ±8 

220 ±9 

LQCD [22] 

278 ± 13 ± 10 

311 ±9 



LQCD [23] 


274 ±6 



LQCD [24] 



175 ±6 

213 ±7 

LQCD (25^ 

245 ± 20 

272 ± 16 

196 ± 24 

229 ± 20 

RPM [30] 

310 

315 

219 

251 

FCM [32] 

273 ± 13 

307 ± 18 

200 ± 10 

230 ± 12 

LFQM [34] 

245^34 

272^ 

196^27 

229^31 

This work 

263 ± 21 

308 ± 21 

213 ± 18 

255 ± 19 


Table 4: The decay constants of the heavy vector mesons from the some theoretical calculations. 



/d* (MeV) 

f DZn { MeV) 

Jb* (MeV) 

/B: n (MeV) 

QCDSR [H] 


128 ± 13 



LQCD [26| 

360 ± 90 

340 ± 110 



This work 

373 ± 19 

333 ± 20 

281 ± 14 

274 ± 13 


Table 5: The decay constants of the heavy scalar mesons from the some theoretical calculations. 


15 































































5m S /p/v/A (MeV) 

<>fs/p/v/A( MeV) 

D (O’) 

+42 

-9 (-11) 

D s (0-) 

+34 

-6 (-8) 

23* (1-) 

+22 

-6 (-8) 

D* s (I") 

+12 

-3 (-5) 

D* 0 (0+) 

-43 

+2 (+6) 

D*s o (0 + ) 

-44 

+1 (+6) 

Di (1+) 

-5 

+1 (+3) 

D s i (1+) 

-2 

+1 (+2) 

23 (O') 

+2 

-3 (-4) 

B s (0-) 

+1 

-2 (-3) 

23* (1") 

+0 

-3 (-3) 

B* (I”) 

+0 

-2 (-2) 

23 0 * (0+) 

-2 

+1 (+3) 

B* s o (0+) 

-1 

+1 (+3) 

Bi (1 + ) 

-0 

+3 (+3) 

B sl (1+) 

-0 

+2 (+2) 


Table 6: The shifts of the masses and decay constants of the heavy-light mesons when the 
perturbative 0(a s ) corrections to the quark condensates are turned off. We can re-choose the Borel 
windows to warrant the mass-shifts Smg/p/v/A = 0, the resulting shifts of the decay constants are 
shown in the bracket. The +0 (—0) denotes the value 0 < 5m < 1 MeV (—1 MeV < 5m < 0). 


we study the masses and decay constants of the pseudoscalar, scalar, vector and axial-vector heavy- 
light mesons with the QCD sum rules in a systematic way. In calculations, we take the MS masses 
and take into account the perturbative 0(a s ) corrections. The masses of the observed heavy-light 
mesons (23,23*), (. D S ,D* S ), (23^(2400), 2?r(2430)), (£>; 0 (2317),£> a i(2460)), (23, 23*), ( B Sl B* s ) can 
be well reproduced, while the predictions for the masses of the (23g,23i) and (23* 0 , 23 s i) can be 
confronted with the experimental data in the futures. Up to the order 0(a s ), the QCD sum rules 
cannot lead to satisfactory values for the fp, fp a and fojfa compared to the experimental data. 
We have to take into account the perturbative 0(aj.) corrections by choosing the pole masses, then 
the experimental data can be well reproduced. The present predictions for the decay constants of 
the heavy-light pseudoscalar, scalar, vector and axial-vector mesons have many phenomenological 
applications in studying the semi-leptonic and leptonic decays of the heavy-light mesons. 
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